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We review our recent studies on selective solvation effects in phase separation in polar binary 
mixtures with a small amount of solutes. Such hydrophilic or hydrophobic particles are preferentially 
attracted to one of the solvent components. We discuss the role of antagonistic salt composed 
of hydrophilic and hydrophobic ions, which undergo microphase separation at water-oil interfaces 
leading to mesophases. We then discuss phase separation induced by a strong selective solvent 
above a critical solute density n p , which occurs far from the solvent coexistence curve. We also give 
theories of ionic surfactant systems and weakly ionized polyelectrolytes including solvation among 
charged particles and polar molecules. We point out that the Gibbs formula for the surface tension 
needs to include an electrostatic contribution in the presence of an electric double layer. 



I. INTRODUCTION 

In soft matter physics, much attention has been paid 
to the consequences of the Coulombic interaction among 
charged objects, such as small ions, charged colloids, 
charged gels, and polyelectrolytesl"^. However, not 
enough effort has been made on solvation effects among 
solutes (including hydrophobic particles) and polar sol- 
vent molecules-^"^. Solvation is also called hydration 
for water and for aqueous mixtures. In mixtures of a 
water-like fluid and a less polar fluid (including poly- 
mer solutions), the solvation is preferential or selective, 
depending on whether the solute is hydrophilic or hy- 
drophobic. See Fig.l for its illustration. The typical 
solvation free energy much exceeds the thermal energy 
krjT per solute particle. Hence selective solvation should 
strongly influence phase behavior or even induce a new 
phase transition. In experiments on aqueous mixtures, 
it is well known that a small amount of salt drastically 
alters phase behavioiEE^. In biology, preferential inter- 
actions between w ater a nd cosolvents with proteins are of 
crutial import anc c! 18 * 19 l Thus selective solvation is rele- 
vant in diverse fields, but its understanding from physics 
is still in its infancy. 

Around 1980, Nabutovskii et a/! 20 | 21 l proposed a pos- 
sibility of mesophases in electrolytes from a coupling 
between the composition and the charge density in 
the free energy. In aqueous mixtures, s uch a cou- 
pling originates from the selective solvatio n 122 * 23 !. It is 
in many cases very strong, as suggested by data of 
the Gibbs transfer free energy in electrochemistry (see 
Sec. 2). Recently, several theoretical groups have pro- 
posed Ginzburg-Landau theories on the solvation in mix- 
ture solvents for electrolytes^ 2 ^ 221 , polyelectrolytes^ 21 ^, 
and ionic surfactants^-. In soft matter physics, such 
coarse-grained approaches have been used to un derstand 
cooperative effects on mesoscopic scales^^l, though 
they are inaccurate on the angstrom scale. They are even 
more usuful when selective solvation comes into play in 
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(a) hydration in water (b) hydration in water-nitrobenzene 



FIG. 1: Illustration of hydration of Na surrounded by a 
shell composed of water molecules in (a) pure water and (b) 
water-nitrobenzene. The solvation chemical potential of Na + 
is higher for (b) than for (a). 



the strong coupling limit. This review presents such ex- 
amples found in our recent research. 

1.1 Antagonistic salt. An antagonistic salt con- 
sists of hydrophilic and hydrophobic ions. An example 
is sodium tetraphenylborate NaBPhj, which dissociates 
into hydrophilic Na + and hydrophobic BPh^ . The lat- 
ter ion consists of four phenyl rings bonded to an ion- 
ized boron. Such ion pairs in aqueous mixtures behave 
antagonistically in the presence of composition hetero- 
geneity, (i) Around a water-oil interface, they undergo 
microphase separation on the scale of the Debye screen- 
ing length while homogeneity holds far from the 
interface to satisfy the charge neutrality (see the right 
bottom plate in Fig. 2). This unique ion distribution pro- 
duces a large elec tric double layer and a large Galvani 
potential differenaP^SESI We found that this ion dis- 
tribution serves to much decrea se the surface tensiorP^, 
in agreement with experiments^^. From x-ray re- 
flectivity measurements, Luo et alW^ determined such 
ion distributions around a water-nitrobenzene(NB) in- 
terface by adding BPh^ and two species of hydrophilic 
ions, (ii) In the vicinity of the solvent criticality, antag- 
onistic ion pairs interact differently with water-rich and 
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oil-rich composition fluctuations, leading to mesophases 
(charge density wave s). In accord with this prediction, 
Sadakane et aZpUMl added a small amount of NaBPh 4 
to a near-critical mixture of D2O and 3-methylpyridine 
(3MP) to find a peak at an intermediate wave number 
g m (^ 0.1 A -1 ~ k) in the intensity of small-angle neu- 
tron scattering. The peak height was much enhanced 
with formation of periodic structures, (iii) Moreover, 
Sadakane et al. observed multi-lamellar (onion) struc- 
tures at small volume fractions of 3MP (in D2 0-rich sol- 
vent) far from the criticality^il, where BPhJ and solvat- 
ing 3MP form charged lamellae. These findings demon- 
strate very strong hydrophobicity of BPI14 . (iv) An- 
other interesting phenomenon is spontaneous emulsifica- 
tion (fo rmati on of small water droplets) at a water-NB 
interface 1 42 * 43 ^. It was observed when a large pure water 
droplet was pushed into a cell containing NB and antag- 
onistic salt (tetraalkylammonium chloride). This insta- 
bility was caused by ion transport through the interface. 

1.2 Precipitation due to selective solvation. 
Many experimental groups have detected large-scale, 
long-lived heterogeneities (aggregates or domains) emerg- 
ing with addition of a hydrophilic salt or a hydropho- 
bic solute in one-phase states of aqueous mixturea^nSS. 
Their small diffusion constants indicate that their typical 
size is of order 10 3 A at very small volume fractions. In 
two-phase states, they also observed a third phase visi- 
ble as a thin solid-like plate at a liquid-liquid interface 
in two-phase states^. In our recent theory^, for suffi- 
ciently strong solvation preference, a selective solute can 
induce formation of domains rich in the selected com- 
ponent even very far from the solvent coexistence curve. 
This phenomenon occurs when the volume fraction of the 
selected component is relatively small. If it is a majority 
component, its aggregation is not needed. This precipita- 
tion phenomenon should be widely observable for various 
combinations of solutes and mixture solvents. 

1.3 Selective hydrogen bonding. Hydrogen bond- 
ing is of primary importance in the phase behavior of 
soft matter. In particular, using statistical-mechanical 
theories, the origin of closed-loop coexistence curve s was 
ascribed to the hydrogen bonding for liquid mixtures ^ 2 * 53 ! 
and for polymer solution ^ 54 * 55 !. Interestingly, water it- 
self can be a selective solute triggering phase separa- 
tion when the hydrogen bonding differs significantly be- 
tween the two compo nents , as observed in a mixture of 
methanol-cyclohexaneSSEZI. More drastically, even water 
absorbed from air changed the phase behavior in films of 
polystyrene(PS)- polyvinylmethylether(PVME^H. That 
is, a small amount of water induces precipitation of 
PVME-rich domains. For block copolymers, similar pre- 
cipitation of micelles can well be expected when a small 
amount of water is added. 

1.4 Ionic surfactant. Surfactant molecules are 
strongly trapped at an interface due to the a mph iphilic 
interaction even if their bulk density is very low^S. They 
can thus efficiently reduce the surface tension, giving rise 
to various mesoscopic structures. However, most theo- 



retical studies have treated nonionic surfactants, while 
ionic surfactants are important in biology and technol- 
ogy. In this review, we also discuss selective solvation in 
systems of ionic surfactants, counterions, and added ions 
in water-oipSl. We shall see that the adsorption behavior 
strongly depends on the selective solvation. 

1.5 Polyelectrolytes. Polyelectrolytes are already 
very complex because of the electrostatic interaction 
among charged particles (ionized monomers and mobile 
ions)pffil Furthermore, we should take into account two 
ingredients^, which have not yet attracted enough atten- 
tion. First, the dissociation (or ionization) on the chains 
should be treated as a chemical reaction in many poly- 
electrolytes containing weak acidic monomera^HSU Then 
the degree of ionization is a space-dependent annealed 
variable. Second, the solvation effects should come into 
play because the solvent molecules and the charged par- 
ticles interact via ion-dipole interaction. Many polymers 
themselves are hydrophobic and become hydrophilic with 
progress of ionization in water. This is because the de- 
crease of the free energy upon ionization is very large. 
It is also worth noting that the selective solvation ef- 
fect can be dramatic in mixture solvents^. As an exam- 
ple, precipitation of DNA has been observed with addi- 
tion of ethanol in water^HsI] where the ethanol added 
is excluded from condensed DNA, suggesting solvation- 
induced wetting of DNA by water. In polyelectrolyte 
solutions, macroscopic phase separation and mesophase 
formation can both take place, sensitively depending on 
many parameters. 

The organization of this paper is as follows. In Sec. 2, 
we will present the background of the solvation on the 
basis of some experiments. In Sec. 3, we will explain a 
Ginzburg-Landau model for electrolytes accounting for 
selective solvation. In Sec. 4, we will treat ionic surfac- 
tants by introducing the amphiphilic interaction together 
with the solvation interaction. In Sec. 5, we will examine 
precipitation induced by a strong selective solute, where a 
simulation of the precipitation dynamics will also be pre- 
sented. In Sec. 6, we will give a Ginzburg-Landau model 
for weakly ionized polyelectrolytes accounting for the ion- 
ization fluctuations and the solvation interaction. In Ap- 
pendix A, we will give a statistical theory of hydrophilic 
solvation at small water contents in oil. 



II. BACKGROUND OF SELECTIVE 
SOLVATION OF IONS 

2.1 Hydrophilic ions in aqueous mixtures. Sev- 
eral water molecules form a solvation shell surrounding 
a small ion via ion-dipole interaction^, as in Fig.l. Clus- 
ter structures produced by sol vation have been observed 
by mass spectrometric analysis^^. Here we mention an 
experiment by Osakai et aW^, which demonstrated the 
presence of a solvation shell in a water-NB mixture in 
two-phase coexistence. They measured the amount of 
water extracted together with hydrophilic ions in a NB- 



3 



rich region coexisting with a salted water-rich region. A 
water content of 0.168M (the water volume fraction <j> 
being 0.003) was already present in the NB-rich phase 
without ions. They estimated the average number of sol- 
vating water molecules in the NB-rich phase to be 4 for 
Na + , 6 for Li + , and 15 for Ca 2+ per ion. Thus, when 
a hydrophilic ion moves from a water-rich region to an 
oil-rich region across an interface, a considerable frac- 
tion of water molecules solvating it remain attached to 
it. Furthermore, using proton NMR spectroscopy. Osakai 
et alW studied successive formation of complex struc- 
tures of anions X~ (such as Cl _ and Br~) and water 
molecules by gradually increasing the water content in 
NB. This hydration reaction is schematically written as 
X"(H 2 0) m _ a + H 2 <— ► X-(H 2 0) m (m = 1,2,3,...). 
For Br~ , these clusters are appreciable for water content 
larger than 0.1M or for water volume fraction <p exceeding 
0.002. 

In Appendix A, we will calculate the statistical distri- 
bution of clusters composed of ions and polar molecules. 
Let the free energy typically decrease by upon bind- 
ing of a polar molecule to a hydrophilic ion of species i. 
A well-defined solvation shell is formed for 3> fcs? 1 , 
where the water volume fraction 4> needs to satisfy 

4> > 4>i ol ~ exp(-e bi /fc B T). (2.1) 

For <fi -c 4>l oi there is almost no solvation. The crossover 
volume fraction </>* ol is very small for strongly hydrophilic 
ions with ebi 3> ksT. For Br~ in water-NB, we estimate 
4>l ol ~ 0.002 from the experiment by Osakai et al. as 
discussed above^. 

2.2 Hydrophobic particles. Hydrophobic objects 
are ubiquitous in nature, which repel water because 
of the strong attrac tion among hydrogen-bonded water 
molecules themselved 7 -* 1 ^. Hyd rophobic particles tend to 
form aggregates in watei^^a and are more soluble in oil 
than in water. They can be either neutral or charged. A 
widely used hydrophobic anion is BPh^ . In water, a large 
hydrophobic particle (> lnm) is even in a cavity sepa- 
rating the particle surface from water^. In a water-oil 
mixture, on the other hand, hydrophobic particles should 
be in contact with oil molecules instead. This attraction 
can produce significant composition heterogeneities on 
mesoscopic scales around hydrophobic o bjects , which in- 
deed takes place around protein surfaces^SDSl. 

2.3 Solvation chemical potential. We introduce a 
solvation chemical potential fJ? so i(<fr) in the dilute limit of 
solute species i. It is the solvation part of the chemical 
potential of one particle (see Eq.(Bl) in Appendix B). 
It is a statistical average over the thermal fluctuations 
of the molecular configurations. In mixture solvents, it 
depends on the ambient water volume fraction <f>. For 
planar surfaces or large particles (such as proteins), we 
may consider the solvation free energy per unit area. 

BorrPl calculated the polarization energy of a polar 
fluid around a hydrophilic ion with charge Z^e using con- 
tinuum electrostatics to obtain the classic formula, 

(Msol)Born = -{Z?e? /2Ri)Q. - 1/e). (2.2) 



The contribution without polarization (e = 1) or in vac- 
uum is subtracted and the (f> dependence here arises from 
that of the dielectric constant e = e(0) (see Eq.(3.5) be- 
low). The lower cutoff Ri is called the Born radi us, which 
is on the order of lA for small metallic ion d 7 l 74 l The hy- 
drophilic solvation is stronger for smaller ions, since it 
arises from the ion-dipole interaction. In this original 
formula, neglected are the formation of a solvation shell, 
the density and composition changes (electrostriction) , 
and the nonlinear dielectric effect. 

For mixture solvents, the binding free energy between 
a hydrophilic ion and a polar molecule is estimated from 
the Born formula (2.1) as £bi ~ — <9[(A t soi) B °m]/<9</> op22l 

e bi ~ Zje 2 e 1 /R l e 2 = k B TZ 2 £ B£l /R t e, (2.3) 

where E\ — de/d<fi and Ib — e 2 /kBTe is the Bjerrum 
length (~ 7A for water at room temperatures). For 
e% ~ e, a well-defined shell appears for Z 2 £ B 3> Ri- 
The solvation chemical potential /4 o i(0) of hydrophilic 
ions in water-oil should largely decrease to negative val- 
ues in the narrow range < <j> < 0* ol , as will be de- 
scribed in Appendix A. In the wide composition range 
(2.1), its composition dependence is still very strong such 
that \dnl ol /d(j)\ ~> k B T holds (see the next subsection). 
The solubility of hydrophilic ions should also increase 
abruptly in the narrow range 4> < 0* ol with addition of 
water to oil. Therefore, solubility measurements of hy- 
drophilic ions would be informative at very small water 
contents in oil. 

In sharp contrast, /i* ol of a neutral hydrophobic par- 
ticle increases with increasing the particle radius R in 
wateM It is roughly proportional to the surface area 
AttR 2 for R > lnm and is estimated to be about 100k B T 
at R ~ lnm. In water-oil solvents, on the other hand, it 
is not easy to estimate the <\> dependence of /4 ol . How- 
ever, /i* ol should strongly increase with increasing the wa- 
ter composition <f>, since hydrophobic particles (including 
ions) effectively attract oil molecule d 19 " 1 !. 

2.4 Gibbs transfer free energy. We consider a 
liquid-liquid interface between a polar (water-rich) phase 
a and a less polar (oil-rich) phase ft with bulk composi- 
tions <f> a and <f)p with <p a > (f>p. The solvation chemical 
potential fi l sol ((j)) takes different values in the two phases 
due to its composition dependence. So we define 

^i^vUM-ulM- (2.4) 

In electrochemistry^- 69 ! 75 * ^, the difference of the solvation 
free energies between two coexisting phases is called the 
standard Gibbs transfer free energy denoted by AG' a ^ for 
each ion species i. Since it is usually measured in units 
of kJ per mole, its division by the Avogadro number A^^ 
gives our A^ a/3 or A/x^ = AG^/Na- With a being the 
water-rich phase, A/i^ is positive for hydrophilic ions 
and is negative for hydrophobic ions from its definition 
(2.4). See Appendix B for relations between A/j,^ and 
other interface quantities. 

For ions, most data of AG^ are at present limited on 
water _ N g69iZ61ZU and wat er-l,2-dichloroethane(EDC)P 
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at room temperatures, where the dielectric constant of 
NB (~ 35) is larger than that of EDC (~ 10). They then 
yield the ratio A/i^/^T. In the case of water-NB, it is 
13.6 for Na+, 27.1 for Ca 2+ , and 11.3 for Br~ as exam- 
ples of hydrophilic ions, while it is —14.4 for hydrophobic 
BPhJ. In the case of water-EDC, it is 22, 7 for Na+ and 
17.5 for Br~, while it is -14.1 for BPhJ . The ampli- 
tude \An l a p\/k B T for hydrophilic ions is larger for EDC 
than for NB and is very large for multivalent ions. In- 
terestingly, A/x^o for H + (more precisely hydronium ions 
H3O" 1 ") assumes positive values close to those for Na + in 
these two mixtures. 

In these experiments on ions, a hydration shell should 
have been formed around hydrophilic ions even in 
the water-poor phase j3 (presumably not completeljESj) . 
From Eq.(2.1) <j>p should be exceeding the crossover vol- 
ume fraction 4> l sol . The data of the Gibbs transfer free 
energy quantitatively demonstrate very strong selective 
solvation even in the range (2.1) or after the shell forma- 
tion. 

On the other hand, neutral hydrophobic particles are 
less soluble in a water-rich phase a than in an oil-rich 
phase p. Their chemical potential is given by //j = 
fcsTln(niAf) + Msol(^)) where * represents the particle 
species and A^ is the thermal de Broglie wavelength. 
From homogeneity of /Lt« across an interface, we obtain 
the ratio of their equilibrium bulk densities as 

n0/n ia = exp[-A^/k B T], (2.5) 
where Afj, l a/3 < from the definition (2.4). 



III. GINZBURG-LANDAU THEORY OF 
MIXTURE ELECTROLYTES 

3.1 Electrostatic and solvation interactions. We 

present a Ginzburg-Landau free energy F for a polar bi- 
nary mixture (water-oil) containing a small amount of 
a monovalent salt [Z\ = 1, Z2 = —1). The multivalent 
case should be studied separately. The ions are dilute and 
their volume fractions are negligibl e, so we are allowed to 
neglect the formation of ion clusters ^ 79 * 80 !. The variables 
cj>, m, and ri2 are coarse-grained ones varying smoothly 
on the molecular sca le. F or simplicity, we also neglect 
the image interactio n 81 ^ 82 ! At a water-air interface the 
image interaction serves to push ions into the water re- 
gion. However, hydrophilic ions are already strongly de- 
pleted from an interface due to their position-dependent 
hydration^]. See our previous analysis^ for relative im- 
portance between the image interaction and the solvation 
interactiion at a liquid-liquid interface. 

The free energy F is the space integral of the free en- 
ergy density of the form, 



The first term / to t depends on 0, ni, and ri2 as 



F = dr 



fu 



1 



C|V0| s 



8tt 



(3.1) 



n, 



ln(n;A?) -l-gi<f> 



(3.2) 



In this paper, the solvent molecular volumes of the two 
components are assumed to take a common value Vq, 
though they are often very different in real bina ry mix- 
tures. Then / is of the Bragg- Williams forr 



vpf 
knT 



= c)>hx(t> + (1 - (£) ln(l - 0) + X cb(l - 0), (3.3) 



where \ is the interaction parameter dependent on T. 
The critical value of x is 2 without ions. The Xi = 
h(2n/m i k B T) 1 ^ 2 in Eq.(3.2) is the thermal de Broglie 
wavelength of the species i with m t being the molecular 
mass. The g\ and 32 are the solvation coupling con- 
stants. In addition, the coefficient C in the gradient part 
of Eq.(3.1) remains an arbitrary constant. To explain ex- 
periments, however, it is desirable to determine C from 
the surface tension data or from the scattering data. 

In the electrostatic part of Eq.(3.1), the electric field is 
written as E = — V<i>. The electric potential $ satisfies 
the Poisson equation, 



V • eV$ = 471-p. 



(3.4) 



The dielectric constant e is assumed to depend on <f> as 



e{(f>) = £ + en 



(3.5) 



where £0 and £1 are positive constants. Though there 
is no reliable theory of e(</>) for a polar mixture, a lin- 
ear composition dependence of e(4>) was observed by De- 
bye and Kleboth for a mixture of nitrobenzene-2,2,4- 
trimethylpentane— . In addition, the form of the elec- 
trostatic part of the free energy density depends on the 
experimental method^. Our form in Eq.(3.1) follows if 
we insert the fluid between parallel plates and fix the 
charge densities on the two plate surfaces. 

We explain the solvation terms in f to t in more detail. 
They follow if Msoi(^) (* = 1' 2) depend on cj> linearly as 



MsolW =Ai~ k B Tgi4>. 



(3.6) 



Here the first term Aj is a constant yielding a contribu- 
tion linear with respect to n, in / to t, so it is irrelevant 
at constant ion numbers. The second term gives rise to 
the solvation coupling in f to t- In this approximation, 
gi > for hydrophilic ions and gi < for hydrophobic 
ions. The difference of the solvation chemical potentials 
in two-phase coexistence in Eq.(2.4) is given by 



A/4 ol = k B T 9i A<j>, 



(3.7) 



where A<fr — <f> a — 4>p is the composition difference. From 
Eqs.(B4) and (B5), the Galvani potential difference is 



A$ = k B T{ 9l - <7 2 )A0/2e, 



(3.8) 
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and the ion reduction factor is 

nxp/ni a = n 2 p/n 2a = exp[-(#i + ,g 2 )A0/2]. (3.9) 

The discussion in subsection 2.3 indicates gi ~ 14 (23) 
for Na + ions and gt — -14 (-14) for BPhJ in water-NB 
(water-EDC) at 300K. For multivalent ions gi can be very 
large (pj ~ 27 for Ca 2+ in water-NB). The linear form 
(3.6) is adopted for the mathematical simplicity and is 
valid for <j> > <fi z sol after the solvation shell formation (see 
Eq.(2.1)). The results in Appendix A suggest a more 
complicated functional form of fJ? so i(4>)- 

In equilibrium, we require the homogeneity of the 
chemical potentials h — 5F/54> and /ij = SF/Srii. Here, 



/' - cv 2 



^E 2 
8ir 



k B T^2gini, (3.10) 



Hi = fc B T[ln(n l A^) - grf] + Zie<&, 



(3.11) 



where /' = df/d<f>. The ion distributions are expressed 
in terms of <f) and $ in the modified Poisson-Boltzmann 
relations^, 



Hi = n°exp[gi4> - Zie®/k B T]. 



(3.12) 



The coefficients arc determined from the conservation 
of the ion numbers, (rii) = V^ 1 J drriiir) — uq, where 
(• • •) = V^ 1 J dr(- ■ •) denotes the space average with V 
being the cell volume. The average uq = (m) = (712) is 
a given constant density in the monovalent case. 

It is worth noting that a similar Ginzburg-Landau free 
energy was proposed by Aerov et al?^* for mixtures of 
ionic and nonionic liquids composed of anions, cations, 
and water-like molecules. In such mixtures, the inter- 
actions among neutral molecules and ions can be pref- 
erential, leading to mesophase formation, as has been 
predicted also by molecular dynamic simulations^. 

3.2 Structure factors and effective ion-ion inter- 
action in one-phase states. The simplest application 
of our model is to calculate the structure factors of the 
the composition and the ion densities in one-phase states. 
They can be measured by scattering experiments. 

We superimpose small deviations 8<f){r) = (f>(r) — {<p) 
and 5rii(r) = Uj(r) — n on the averages (</>) and no- 
The monovalent case {Z\ = 1,Z 2 = —1) is treated. As 
thermal fluctuations, the statistical distributions of 8<j>(r) 
and Srii are Gaussian in the mean field theory. We may 
neglect the composition-dependence of e for such small 
deviations. We calculate the following, 

S(q) = (\(t>q\ 2 ) c , G l3 {q) = (n t qn* q )Jn , 

C(q) = {\pq\ 2 )Je 2 n , (3.13) 

where <f>q, mq, and pq are the Fourier components of 6<j>, 
Srii, and the charge density p = e(n\ — n 2 ) with wave 
vector q and (• • -) denotes taking the thermal average. 
We introduce the Bjerrum length l B = e 2 /ek B T and the 
Debye wave number k = (Sirisno) 1 ^ 2 . 



First, the inverse of S(q) is written a d 22 * 23 * 

2 



1 - / , ^ 2 n Cq 



1 



1> 2 



, (3.14) 



where f = f"/k B T with /" = d 2 f/d<j) 2 . The second 
term is large for large (171 + g 2 ) 2 even for small average 
ion density no, giving rise to a large shift of the spin- 
odal curve. If g\ ~ g 2 ~ 15, this factor is of order 10 3 . 
In the previous experiments^^, the shift of the coexis- 
tence curve is typically a few Kelvins with addition of a 
10~ 3 mole fraction of a hydrophilic salt like NaCl. The 
parameter 7 P in the third term represents asymmetry of 
the solvation of the two ion species and is defined by 



7 P = {k B T/lQ-KCt B ) 1 l 2 \g 1 -g 2 



(3.15) 



If the right hand side of Eq.(3.14) is expanded with re- 
spect to q 2 , the coefficient of q 2 is C(l — ~/ 2 )/k B T. Thus 
a Lifshitz point is realized at 7 P = 1. For 7 P > 1, S(q) 
has a peak at an intermediate wave number, 



<7m = (7p - 1) 



1/2, 



(3.16) 



The peak height S(q m ) and the long wavelength limit 
5(0) are related by 



1/S(q m ) - 1/5(0) - C( 7p - iy K 2 /k B T. 



(3.17) 



A mesophase appears with decreasing f or increasing x, 
as observed by Sadakane et alW^. In our mean-field the- 
ory, the criticality of a binary mixture disappears if a 
salt with 7 P > 1 is added however small its content 
is. Very close to the solvent criticality, Sadakane et 
alW^. recently measured anomalous scattering from D 2 0- 
3MP-NaBPh 4 stronger than that from D 2 0-3MP with- 
out NaBPh.4. There, the observed scattering amplitude 
is not well described by our S(q) in Eq.(3.14), requiring 
more improvement. 

Second, retaining the fluctuations of the ion densities, 
we eliminate the composition fluctuations in F to obtain 
the effective interactions among the ions mediated by the 
composition fluctuations. The resultant free energy of 
ions is written as^ 

■Fion = / dr k B T ni ln(^A 3 ) 

J i 

+i fdrfdr'Y f V ij (\r-r'\)5n i (r)Sn j (r'). (3.18) 

id 

The effective interaction potentials V«(r) are given by 

V ij {r)=Z i Z s --^-e- r ^, (3.19) 
er A r 

where Z\ and Z 2 are ±1 in the monovalent case, A = 
4nC/(k B T) 2 , and £ = (C/f) 1/2 is the correlation length. 
The second term in Eq.(3.19) arises from the selective 
solvation and is effective in the range a < r < £ and 
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can be increasingly important on approaching the solvent 
criticality (for £ ^> a). It is attractive among the ions of 
the same species (i — j) dominating over the Coulomb 
repulsion for 



g\ > 4irC£ B /k B T. 



(3.20) 



Under the above condition there should be a tendency of 
ion aggregation of the same species. In the antagonistic 
case (<?i<72 < 0), the cations and anions can repel one 
another in the range a < r < £ for 



\g l92 \>^Cl B /k B T, 



(3.21) 



under which charge density waves are triggered near the 
solvent criticality. 

The ionic structure factors can readily be calculated 
from Eqs.(3.18) and (3.19). Some calculations give^l 



Gu(q) 
Go(q) 



1 



GMq) = g 



n S(q) 
n 



9l 



(.gi - g 2 ) 2 2g 2 u 



2{u+l) 2m + 1 



S(q)(gi + .92) 



4 C(q), 



C(q) = 



2u 



al d9l ~.92) 2 2 

n oS{q)-, — r^-u , 

(u + iy 



where u = q 2 /k 2 and 



Go(q) = 



u+1/2 _ q 2 + K 2 /2 
u + 1 g 2 + k 2 



(3.22) 



(3.23) 



is the structure factor for the cations (or for the anions) 
divided by no in the absence of solvation. The solva- 
tion parts in Eq.(3.22) are all proportional to noS(q), 
where S(q) is given by Eq.(3.14). The Coulomb inter- 
action suppresses large-scale charge-density fluctuations, 
so C(q) tends to zero as q — > 0. 

It should be noted that de Gennes^ derived the 
effective interaction among monomers on a chain (oc 
— e~ r l^ jr) mediated by the composition fluctuations in 
a mixture solvent. He then predicted anomalous size be- 
havior of a chain near the solvent criticality. However, 
as shown in Eq.(3.19), the effective interaction is much 
more amplified among charged particles than among neu- 
tral particles. This indicates importance of the selective 
solvation for a charged polymer in a mixture solvent, even 
leading to a prewetting transition around a chaiipH We 
also point out that an attractive interaction arises among 
charged colloid particles due to the selective solvation in 
a mixture solvent, on which we will report shortly. 

3.3 Liquid-liquid interface profiles. The second 
application is to calculate a one-dimensional liquid-liquid 
interface at z — zq taking the z axis in its normal direc- 
tion, where (j) — > (f> a in water-rich phase a (z — Zq — > —00) 
and <j> — > <j)fj in oil- rich phase f3 (z — zq — > 00). In 
Fig. 2, we give numerical results of typical interface pro- 



files, where we measure space in units of a 



- 1/3 



and set 



C = 2.5k B T/a 2 , e 2 /s k B T = 3a, and e\ = e . In these 
examples, the correlation length £ is shorter than the 



Debye lengths K a x and /to 1 in the two phases. However, 

near the solvent criticality, £ grows above k" 1 and n^ 1 
and we encounter another regime, which is not treated 
in this review. 
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X = 2 2 
/ aC = 2.5 

gl = 1 - 
g 2 - 13 
lO^orti. = 6.27 
e 2 /E B k B T = 3a - 






£1 = E(j 

/ 




*(z)/|A*| 









1 


X = 2.2 




aC- 2.5 




91 =7 




92 = 13 




llrWtio = 6.27 " 




e'/toksT -3a 




ei — E a 











400 

zja 



4C0 

zja 



-*(js)/|A*| 



0.8 - fa) 



X = 2.45 
aC - 2.5 
- 9i ~ -ff2 = 15 
lOWli. = 1-63 

e'/a,k B T = 3o 
ei =so 



220 



240 260 

zja 



X = 2.45 




uC - 2.5 




- 3i = -32 = 15 . 




10'tionia = 1.63 | 




e'/EakuT = 3a 




£l -£c 




«0"2(z) J 





240 2S0 

zja 



FIG. 2: Top: normalized potential $(z)/A$ (left), compo- 
sition <f>(z) (left), and normalized ion densities voni(z) and 
vori2(z) (right on a semilogarithmic scale) for hydrophilic ion 
pairs, where \ = 2.2, gi = 7, gi = 13, Vo7ii Q = vonip = 
4 x 10~ 4 , and eA<3> = — 1.55fcsT. Bottom: those for an- 
tagonistic ion pairs, where \ = 2.45, gi = 15, $2 = —15, 
1.63 x 10" 4 , and eA$ = 10.3fc s T. In 



v n lc 



vamp 



these cases e jeoksT = 3a and ei 



The upper plates give <p(z), $(z), n 1 (z), and 712(2) 
for hydrophilic ion pairs with gi = 7 and g 2 = 13 at 



X = 2.2 and n\ a — n 2a = 4 x 10 



The ion reduc- 



tion factor in Eq.(3.9) is 0.0057. The potential <& varies 
mostly in the right side (phase j3) on the scale of the De- 
bye length nZ 1 = 67.5a (which is much longer than that 
k^ 1 = 6.1a in phase a). The Galvani potential differ- 
ence A<E> is 0.76k B T/e here. The surface tension here is 
a = 6.53 x 10~ 2 k B T/a 2 and is slightly larger than that 
(j q = 6.22k B T/a 2 without ions (see the next subsection). 

The lower plates display the same quantities for antag- 
onistic ion pairs with 171 = 15 and g 2 = —15 for \ = 2.45 
and iii a = n 2a = 1.67 x IO^^Vq 1 . The anions and the 
cations are undergoing microphase separation at the in- 
terface on the scale of the Debye lengths k~ 1 = 12.3a and 
k~o 1 — 9.70a, resulting in a large electric double layer and 
a large potential drop (~ 10k B T/e). The surface ten- 
sion here is a = 0.0805k B T/a 2 and is about half of that 
(j = 0.155k B T/a 2 without ions. This large decrease in 
a is marked in view of small ri\ a . A large decrease o f the 
surface tension was observed for an antagonistic sail! 37 * 38 ! 

3.4 Surface tension. There have been numerous 
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measurements of the surface tension of an air- water inter- 
face with a salt in the water region. In this case, a lmost 
all salts lead to an increase in the surface tensiorP3SU, 
while acids tend to lower it because hydronium ions are 
trapped at an air- water interface! 88 * 89 ! 

Here, we consider the surface tension of a liquid-liquid 
interface in our Ginzburg-Landau scheme, where ions can 
be present in the two sides of the interface. In equilibrium 
we minimize Q, = J dzoj, where u is the grand potential 
density, 

u = /tot + \c\V4>\ 2 + ^E 2 -h4>-J2 mrH. (3.24) 

i 

Using Eqs.(3.10) and (3.11) we hnd d(cu + p$)/dz = 
where cj>' = d<f>/dz and 4>" = d 2 (f>/dz 2 . Thus, 



C(f) 12 - p$ 



(3.25) 



Since (/>' and p tend to zero far from the interface, uj(z) 
tends to a common constant uj^ as z — > ±oo. T he su rface 
tension a = j dz\ui{z) — lOqo] is then written a d 23 * 24 ! 



dz 



C<j>' 2 



^E 2 
4tt 



2u„ — 2(T , 



(3.26) 



where we introduce the areal densities of the gradient free 
energy and the electrostatic energy as 



dz C(j)' 2 /2, 



dz sE 2 /8tt. 



(3.27) 



The expression a = 2a s is well-known in the Ginzburg- 
Landau theory without the electrostatic interaction^. 



where cto is the surface tension without ions and T is 
the adsorption to the interface. In terms of the total ion 
density n = rii + n 2 , it may be expressed as 



T = dz 



n — n a — ~ M 



(3.29) 



where uk = ti\k + n%K (K = ex, /3), An = n a — n», and 
the integrand tends to zero asz4 ±oo. From Eqs.(3.26) 
and (3.28) <r g is expressed at small ion densities as 



2er„ = a — k B TT + a e 



(3.30) 



In the Gibbs formula (a ^ a - k B TTf^, the electro- 
static contribution — a e is ne glected. However, it is cru- 
cial for antagonistic salt! 24 * 29 * and for ionic surfactant)^. 

In Fig. 3, numerical results of 2er e , 2er g , er, and the com- 
bination (To — k B TT + a c are plotted as functions of the 
bulk ion density n a — ni a + ni a for the antagonistic 
case (?i = —g% = 15, where % = 2.4 and C — 2.5k B T/a 2 . 
The parameter 7 P in Eq.(3.15) exceeds unity (being equal 
to 1.89 for e = 1.5eo). In this example, a g weakly 
depends on n a and is fairly in accord with Eq.(3.30), 
while cr c steeply increases with increasing n a . As a re- 
sult, a e increases up to <7 g , leading to vanishing of a at 
n Q = 8 x lCT 4 Vo 1 . 

We may understand the behavior of a e as a function 
of g\ andg2 by solving the nonlinear Poisson-Boltzmann 
equation^ with an interface at z = 0. That is, away 
from the interface \z\ > £, the normalized potential 
U(z) = e$(z)/k B T obeys 



; U = k% smh(U-U K ) 



(3.31) 



o 
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FIG. 3: Surface quantities a, 2a c , 2a s , and cro — ksTT + a e in 
units of fcsTa -2 vs von a for a monovalent, antagonistic salt 
with gi — —g<i = 15 at X = 2.4. Here 2(T g and oo — ksTT + a e 
are close, supporting Eq.(3.30). Growth of 2a e gives rise to 
vanishing of the surface tension a at von a = 8 x 1CP 4 . 



In our previous worlJMSi!, we obtained the following 
approximate expression for a valid for small ion densities: 



(j — do — k B TT 



in the two phases (K = a, (3), where U a = U{—oo), 
Up = U(oo), and k k — {^Trn K e 2 j 'exksT) 1 / 2 with e K — 
Eq +£\4>k- In solving Eq.(3.31) we assume the continuity 
of the electric induction —ed$/dz at z — (but this 
does not hold in the presence of interfacial orientation of 
molecular dipoles, as will be remarked in the summary 
section). The Poisson-Boltzmann approximation for a e 
is of the form, 



k B T 



2n 

A s (n a /£ Ba )^ 2 



y/l + b 2 + 26cosh(A(7/2) -6-1 

(3.32) 



We should have a c = cr^ B in the thin interface limit 
£ <C nj/ . In the first line, the coefficient b is defined by 



b = (e p/e a y ' 2 exp[-(. 9l + ,g 2 )A0/4], 



(3.33) 



and AU — U a — Up — {g\ — g2)A<fi/2 is the normal- 
ized potential difference calculated from Eq.(3.8). In the 
second line, 4o = e 2 /e a k B T is the Bjerrum length in 
phase a. The second line indicates that the electrostatic 
contribution to the surface tension is negative and is of 



order n^l 2 as n 



away from the solvent critical- 



(3.28) ity, as first predicted by Nicols and Pratt for liquid-liquid 



interfaces^. Remarkably, the surface tension of air- water 
interfaces exhibited the same behavior at very small salt 
densities (known as the Jones-Ray effect pH though it 
has not yet been explained reliably^. 

In the asymptotic limit of antagonistic ion pairs, we 
assume g\ > —52 3> 1, where the coefficient A s in the 
second line of eq.(3.32) grows as 

A s = 7r-V 2 ( e/J / ea )V4exp(|<»|A0/4) (3.34) 

We may also examine the usual case of hydrophilic ion 
pairs in water-oil, where gi and c/2 are both considerably 
larger than unity. In this case A s becomes small as 

A s = (£ /3 /^ C[ ) 1 / 2 [cosh(A{//2) - f] 

xexp[-( ffl + (/2)A0/4]. (3.35) 

1/2 

In this case, the electrostatic contribution — <r c (cx nj ) in 
a could be detected only at extremely small salt densities. 

Analogously, between ionic and nonionic liquids, Aerov 
el a/.^ calculated the surface tension. They showed that 
if the affinities of cations and anions to neutral molecules 
are very different, the surface tension becomes negative. 

3.5 Mesophase formation with antagonistic salt. 
Adding an antagonistic salt with 7 P > I to water-oil, we 
have found instability of one-phase states with increas- 
ing x below Eq.(3.I5) and vanishing of the surface ten- 
sion a with increasing the ion content as in Fig. 3. In 
such thermodynamic instability is induced with 

increasing x at a fixed ion density n = (ni) = (712), 
leading to a mesophase. To examine this ph ase or der- 
ing, we performed two-dimensional simulation ;] 29 ! 30 ! and 
presented an approximate phase diagrarrPH. We here 
present preliminary three-dimensional results. The pat- 
terns to follow re sembl e those in block copolymers and 
surfactant systems^^H. In our case, mesophases emerge 
due to the selective solvation and the Coulomb interac- 
tion without complex molecular structures. Solvation- 
induced mesophase formation can well be expected in 
polyelectrolytes and mixtures of ionic and polar liquids. 

We are interested in slow composition evolution with 
antagonistic ion pairs, so we assume that the ion distri- 
butions are given by the modified Poisson-Boltzman n re- 
lations in Eq.(3.I2). The water composition (f> obeys^l22l 

! + V.(« = ^V*„, (3.36) 

where Lq is the kinetic coefficient and h is defined by 
Eq.(3.I0). Neglecting the acceleration term, we deter- 
mine the velocity field v using the Stokes approximation, 

r] V 2 v = Vp 1 + < />V/i + ^n i V£t i , (3.37) 

i 

where r\§ is the shear viscosity and \x\ are defined by 
Eq.(3.1I), We introduce pi to ensure the incompressibil- 
ity condition V • v = 0. The right hand side of Eq.(3.37) 
is also written as V • II, where II is the stress tensor aris- 
ing from the fluctuations of <j) and rii . Here the total free 




FIG. 4: Composition patterns at t = 5000io for x = 2.1 
and Van = 3 x 10~ 3 with average composition {<f>} being 0.5 
(left) and 0.4 (right) for antagonistic salt with gi — —g% = 10. 
These patterns are nearly pinned. Yellow surfaces are oriented 
to the regions of <j> > 0.5 and blue ones are to those of < 0.5. 




0.0 ^^^w 0.01 



FIG. 5: Cross-sectional profiles of cation n\ (left) and anion 
n-2 (right) in the x-y plane at z — for antagonistic salt using 
data in Fig. 4. The domains are bicontinuous for (<fi) = 0.5 
(top) and droplet-like for {(j>) — 0.4 (bottom). 



energy F in Eq.(3.1) satisfies dF/dt < with these equa- 
tions (if the boundary effect arising from the surface free 
energy is neglected). 

We integrated Eq.(3.36) using the relations (3.4), 
(3.12), and (3.37) on a 64 x 64 x 64 lattice under the 
periodic boundary condition. The system was quenched 
to an unstable state with x = 2.1 at t = 0. Space and 
time are measured in units of a = v^ 3 and to = a 5 /L , 
respectively. Without ions, the diffusion constant of the 
composition is given by Lof" /ksT in one-phase states 
in the long wavelength limit (see Eq.(3.14)). We set 
91 = -g 2 - 10, n - 3 x 10- 3 «o X . C = k B T/2a 2 , 
e 2 /eoksT = 3a, e\ = 0, and ri^L^/kBT = 0.16a 4 . 

In Fig. 4, we show the simulated domain patterns at 
t = 5000to, where we can see a bicontinuous structure for 
(4>) = 0.5 and a droplet structure for (4>) = 0.4. There 
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FIG. 6: Normalized structure factor S(q)/vo for antagonistic 
salt from domain structures in pinned states at t = 5000to, 
where the patterns are bicontinuous at for (</>) = 0.5. A peak 
height S(q m )/vo of each curve much exceed unity. The aver- 
age ion density is vono — 0.001, 0.002, and 0.003 for the three 
curves (from above for q < q m )- 



is almost no further time evolution from this stage. In 
Fig. 5, the ion distributions are displayed for these two 
cases in the x-y plane at z = 0. For (0) = 0.5, the ion 
distributions are peaked at the interfaces forming elec- 
tric double layers (as in the right bottom plate of Fig. 2). 
For ((f)) = 0.4, the anions are broadly distributed in the 
percolated oil region, but we expect formation of electric 
double layers with increasing the domain size also in the 
off-critical condition. In Fig. 6, the structure factor S(q) 
in steady states are plotted for vqUq = 0.001,0.002, and 
0.003. The peak position decreases wi th inc reasing no in 
accord with Eq.(3.16). Sadakane et a/! 39 l 41 l observed the 
structure factor similar to those in Fig. 6. 

Finally, we remark that the thermal noise, which is 
absent in our simulation, should be crucial near the crit- 
icality of low-molecilar-weight solvents. It is needed to 
explain anomalously enhanced composition fluctuations 
induced by NaBPli4 near the solvent criticalitySSl. 



IV. IONIC SURFACTANT WITH 
AMPHIPHILIC AND SOLVATION 
INTERACTIONS 

4.1 Ginzburg-Landau theory. In this section, we 
will give a diffuse-interface model of ionic surfactant^', 
where surfactant molecules are treated as ionized rods. 
Their two ends can stay in very different environments 
(water and oil) if they are longer than the interface thick- 
ness £. In our model, the adsorption of ionic surfac- 
tant molecules and counterions to an oil-water interface 
strongly depends on the selective solvation parameters <?i 
and gi and that the surface tension contains the electro- 
static contribution as in Eqs.(3.26) and (3.28). 

We add a small amount of cationic surfactant, anionic 



counterions in water-oil in the monovalent case. The den- 
sities of water, oil, surfactant, and counterion are ua, 
tib, nx, an d n 2 , respectively. The volume fractions of 
the first three components are 4>a — vo n A, 4>b = vqTIb, 
and Virii, where vq is the common molecular volume of 
water and oil and v\ is the surfactant molecular volume. 
The volume ratio N\ = V\/vq can be large, so we do not 
neglect the surfactant volume fraction, while we neglect 
the counterion volume fraction supposing a small size of 
the counterions. We assume the space-filling condition, 



t>A + 4>B + Villi = 1. 



(4.1) 



Let 2ip = 4>a — 4>b be the composition difference between 
water and oil; then, 



b A = (1 - «mi)/2 + ip, 

b B = (l - «mi)/2 - if). 



(4.2) 



The total free energy F is again expressed as in 
Eq.(3.1). Similarly to Eq.(3.2), the first part reads 



1 



[4>a In (/>a + 4>B In (j>B + X^A^i 



k B T v ' 

+ 2jni[ln(niAf) - gifiiip] ~ n x \nZ a , (4.3) 

i 

The coefficients g\ and g 2 are the solvation parameters 
of the ionic surfactant and the counterions, respectively. 
Though a surfactant molecule is amphiphilic, it can have 
preference to water or oil on the average. The last term 
represents the amphiphilic interaction between the sur- 
factant and the composition. That is, Z a is the partition 
function of a rod-like dipole with its center at the posi- 
tion r. We assume that the surfactant molecules take a 
rod-like shape with a length 2£ considerably longer than 

1/3 

a = v . It is given by the following integral on the 
surface of a sphere with radius I, 



-in 



exp 



Wa,ip(r — £u) — w a tp(r + £u) 



(4.4) 



where u is the unit vector along the rod direction and 
/ d£l represents the integration over the angles of u. The 
two ends of the rod are at r + £u and r — lu under the 
influence of the solvation potentials given by kBTw li tjj( r + 
lu) and —k B Tw a ip(r—£u). The parameter u> a represents 
the strength of the amphiphilic interaction. 

Adsorption is strong for large w a A-0 3> 1, where Aip = 
ip a — ipp(— <pAa — i'Ap) is the difference of tp between the 
two phases a and (3. In the one-dimensional case, all the 
quantities vary along the z axis and Z a is rewritten as 



Z a (z) 



2£ 



exp[w a -0(z - C) - w a ip(z + ()]. (4.5) 



where £ = £u z . In the thin interface limit £ <C £, we place 
the interface at z — to find Z a — 1 for \z\ > £, while 

Z a (z) = 1 + (1 - |z|/£)[cosh( Wa AV0 - 1], (4.6) 
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for \z\ < £. Furthermore, in the dilute limit villi <C 1 and 
without the electrostatic interaction, we have n\(z) — 
rii Q Z a (z) for z < and ni(z) = n±pZ & (z) for z > 0, 
where n\ a and nip — e~ 9lA ^'ni a are the bulk surfactant 
densities. The surfactant adsorption then grows as 



Ti = 



dz[ni(z) — n a ] + / dz[ni(z) — np] 



(n la + nip)£{cosh(w a Aip) - l]/2. 



(4.7) 



However, the steric effect comes into play at the interface 
with increasing the surfactant volume fraction at the in- 
terface (~ TiVi/£). 
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FIG. 7: Profiles for mixtures with cationic surfactant and 
anionic counterions with v\ — 5vo and w a = 12. Top: Wini 
(bold line), <j>A, and <^g- Middle: Voni and uo^2- Bottom: 
e^/ksT exhibiting a maximum at the interface. Here g\ = 4, 
<?2 = 10, and vonia — 10~ 3 (left), while gi = —g% — 8 and 
vonic = 3.6 x 1CP 4 (right). The counterion distribution has 
a peak in the phase a (left) or f3 (right) depending on g-2- 
[From: A. Onuki, Europhys. Lett. 82, 58002 (2008)]. 



4.2 Interface profiles of compositions, ion den- 
sities, and potential. We give typical one-dimensional 
interface profiles varying along the z-axis in Fig. 7. We set 
t'l = 5«o, C = 3ksT/a, x = 3, and e 2 /a£ c /csT = 16/tt. 
The dielectric constant is assumed to be of the form 
e = e c (l + 0.8?/>), where e c is the critical value. Then 




10 v n 1( 



FIG. 8: aa 2 /k B T, (a + a c )a 2 /k B T, and Fa 2 as functions 
of vonia with vi = 5vo, x = 3, and iw a = 12. The curves 
change on a scale of 10 -3 for hydrophilic ion pair g± — 4 and 
32 = 10 (top) and on a scale of 10 -4 for antagonistic ion pair 
gi = —(72 = 8 (bottom). [From: A. Onuki, Europhys. Lett. 
82, 58002 (2008)]. 



e a = 2Ep at x — 3. This figure was produced in the 
presence of the image interaction in our previous workP^l 
(though it is not essential here). 

In Fig. 7, we show the volume fractions <f>A, 4>b, and 
vini = 1 — cf>A — 4>b (top), the ion densities n\ and 
n 2 (middle), and the potential e^/fc^T with $ Q = 
(bottom). In the left, the counterions are more hy- 
drophilic than the cationic surfactant, where g\ = 4 and 
10 leading to V = 0.124a" 2 and a = QM7k B Ta~ 2 



-3„,-l 



In the right plates, the surfac- 



es 

at n la = 10 

tant cations are hydrophilic and the counterions are hy- 
drophobic, where g\ = —g 2 = 8 leading to T = 0.155a -2 
and a = 0.l59k B T a - 2 at n la = 3.6 x lO -4 ^ x . The 
distribution of the surfactant n\ is narrower than that 
of the counterions n 2 - This gives rise to a peak of <E> at 
z = z p , at which E(z p ) oc f^_" dz(ni(z) — n 2 {z)) = 0. 

The adsorption strongly depends on the solvation pa- 
rameters <7i and g 2 . It is much more enhanced for antag- 
onistic ion pairs than for hydrophilic ion pairs. 

4.3 Surface tension. The grand potential density 
u) is again given by Eq.(3.24) and tends to a common 
constant ujrx, as z — > ±oo, though its form is more com- 
plicated. The surface tension a = J dz[ui(z) — Wm] is 
rewritten as Eq.(3.26) and is approximated as Eq.(3.28) 
for small The areal electrostatic-energy density a e 
in Eq.(3.27) is again important in the present case. 
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In Fig. 8, we show a, a + er c , and T as functions of 
vonia at w a — 12, where T is defined as in Eq.(3.28) for 
n = ni+n 2 - In the upper plate, the two species of ions are 
both hydrophilic (gi = 4 and <?2 = 10), while in the lower 
plate the surfactant and the counterions are antagonistic 
(<?i = 8 and 172 = —8). For the latter case, a large electric 
double layer is formed at an interface, leading to a large 
r and a large decrease in a even at very small ri\ a . In the 
present case the Gibbs term k B TT is a few times larger 
than a e . Note that the approximate formula (3.28) may 
be derived also in this case, but it is valid only for very 
small ri\ a in Fig. 8. 

V. PHASE SEPARATION DUE TO STRONG 
SELECTIVE SOLVATION 

5.1 Strongly hydrophilic or hydrophobic solute. 

With addition of a strongly selective solute in a binary 
mixture in one-phase states, we predict precipitation of 
domains composed of the preferred component enriched 
with the solute^l. These precipitation phenomena occur 
both for a hydrophilic salt (such as NaCl) and a neutral 
hydrophobic solut c 44 " 50 [ In our scheme, a very large size 
of the selective solvation parameter is essential. In 
Sees. 2 and 3, we have shown that |<7j| can well exceed 10 
both for hydrophilic and hydrophobic solutes. 

With hydrophilic cations and anions, a charge den- 
sity appears only near the interfaces, shifting the surface 
tension slightly. Thus, in the static aspect of precipita- 
tion, the electrostatic interaction is not essential, while 
fusion of precipitated domains should be suppressed by 
the presence of the electric double layers. We will first 
treat a hydrophilic neutral solute as a third component, 
but the following results are applicable also to a neutral 
hydrophobic solute if water and oil are exchanged. In ad- 
dition, in a numerical example in Fig.ll, we will include 
the electrostatic interaction among hydrophilic ions. 

5.2 Conditions of two phase coexistence. Adding 
a small amount of a highly selective solute in water-oil, 
we assume the following free energy density, 

MM = f(^) + k B Tn[ln(n\ 3 ) - 1 - g<f>]. (5.1) 

This is a general model for a dilute solute. For monova- 
lent electrolytes, this form follows from Eq.(3.11) if there 
is no charge density or if we set 

ni = n 2 = 7i/2, 9 = {9i+ Sa)/2, (5.2) 

The first term /(</>) is assumed to be of the Bragg- 
Williams form (3.3). The A is the thermal de Broglie 
length. The solvation term (cx g) arises from the so- 
lute preference of water over oil (or oil over water) . The 
strength g is assumed to much exceed unity^. We fix the 
amounts of the constituent components in the cell with 
a volume V. Then the averages 4> — (<P) — I dr<f>/V and 
n = (n) = J drn/V are given control parameters as well 
as x- 



In two phase coexistence in equilibrium, let the com- 
position and the solute density be (<p a ,n a ) in phase a 
and (<f>p,rip) in phase /3, where 4> a > 4> > <f)p and 
n a > n > rip. We introduce the chemical potentials 
h = dftot/d<j) and /j = df to t/dn. Equation (5.1) yields 

h = f'(<p) - k B Tgn, (5.3) 
p t /k B T = 1n(nX 3 )-g4>, (5.4) 

where /' = df/4>. The system is linearly stable for 
dh/d(j) - (dh/dnf/dfi/dn > or for 

f'(4>) - k B Tg 2 n > 0, (5.5) 

where /" = d 2 f/d<fi 2 . Spinodal decomposition occurs if 
the left hand side of Eq.(5.5) is negative. 
The homogeneity of (i yields 

n = ne^/(e^), (5.6) 

where (• • •) denotes taking the space average. The bulk 
solute densities are uk = fie 9 ^ K / (e 9 ^) for K — a, (3 in 
two-phase coexistence. In our approximation Eq.(5.6) 
holds even in the interface regions. We write the volume 
fraction of phase a as j a . We then have (e 9 ^) = j a e 9 ^ a + 
(1 — 7 Q )e s ^' 3 in Eq.(5.6). In terms of <f> and n, 7 Q is 
expressed as 

7a = (4> - (t>p)/&<t> = {n - rip)/ An. (5.7) 

where A<p = <j) a — (f)p > and An = n a — np > 0. In 
these expressions we neglect the volume of the interface 
regions. Since n a jnp = e 9 ^ >• 1 for gA<f> ^> 1, the 
solute is much more enriched in phase a than in phase j3. 
Eliminating n using Eq.(5.6), we may express the average 
free energy density as 

(/tot) = (/) - k B Tn\n{(e 9 *)} + A x , (5.8) 

where A\ = k B Tn[\n(n\ 3 ) — 1] is a constant at fixed n. 
In terms of <j) a , 4>p, and 7„, Eq.(5.8) is rewritten as 

(/tot) = [7a/(0a) + (1 - 7a)/(^)] 

~k B Tn\i,[ la e 9 ^ + (1 - 7a)e 9 ^] + A x . (5.9) 

The second term (oc ft) is relevant for large g (even for 
small h). Now we should minimize (/tot) — h[y a (f> a + 
(1 — 7a)</>£ — 4>] with respect to <p a , (ftp, and j a at fixed </>, 
where h appears as the Lagrange multiplier. Then we ob- 
tain the equilibrium conditions of two-phase coexistence, 

h = f{4> a ) ~ k B Tgn a = f(<j>p) ~ k B Tgnp, (5.10) 
fiM - f(<j>p) - k B TAn = hAcf,. (5.11) 

These static relations hold even for ion pairs under 
Eq.(5.2). 

5.3 Numerical results of two phase coexistence. 

In Fig. 9, we give numerical results on the phase behavior 
of (j) a and 4>p in the left and n a and j a in the right 
as functions of %. We set g = 11, (j> = 0.35 and n = 
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FIG. 10: Leit:x P {4>,n) vs </> for three values of n at j = 11, 
which nearly coincide with the asymptotic formula (5.22) 
(dotted line) for <j> < 0.35 and converge to the coexistence 
curve for larger <f>. Coexistence region without ions is in the 
upper region (in right green). Right: von p ((j>, \) vs ^ on a 
semi- logarithmic scale at g = 11, It coincides with the asymp- 
totic formula (5.23) (dotted line) for \ < 2 and tends to zero 
at the coexistence composition, 0.249 (top) or 0.204 (bottom). 
[From: R. Okamoto and A. Onuki, Phys. Rev. E 82, 051501 
(2010).] 



FIG. 9: Left: compositions <f> a and <j>p vs \- Right: semi- 
logarithmic plots of volume fraction 7 a and normalized solute 
density n a /n of the water-rich phase a vs x- Here g = 11, 
4> = 0.35, and ft = 6 x lO^ 4 ^ 1 ( to P)> while 9 = 10 - 5 . = 
0.58, and n = 4 x lO" 3 ^ 1 (bottom). Plotted in the left 
also are the coexisting region without solute (in right green) 
and the spinodal curve with solute (on which the left hand 
side of Eq.(5.5) vanishes). See Fig. 11 for simulation of phase 
separation on point (A) in the left bottom panel. [Upper 
plates are from: R. Okamoto and A. Onuki, Phys. Rev. E 
82, 051501 (2010).] 



6 x 10~ 4 Vq 1 in the top plates and g = 10.5, <f> — 0.58, 
and n — 4 x lO -3 ^ 1 in the bottom plates. The solute 
density is much larger in the latter case. Remarkably, a 
precipitation branch appears in the range, 

X P < X < 2. (5.12) 

The volume fraction j a decreases to zero as x approaches 
the lower bound Xp = Xp(^j "•)• Without solute, the mix- 
ture would be in one-phase states for \ < 2. The pre- 
cipitated domains are solute-rich with <f> a = 1, while 4>p 
is slightly larger than </>. In the left upper plate </> a in- 
creases continuously with decreasing while in the left 
lower plate 4> a jumps at x = 1-937 and hysteresis appears 
in the region 1.927 < x < 1-952. We also plot the spin- 
odal curve, f"((f>) — kBTg 2 n = 0, following from Eq.(5.5). 
Outside this curve, homogeneous states are metastable 
and precipitation can proceed via homogeneous nucle- 
ation in the bulk or via heterogeneous nucleation on hy- 
drophilic surfaces of boundary plates or colloids 32 . Inside 
this curve, the system is linearly unstable and precipita- 
tion occurs via spinodal decomposition. This unstable 
region is expanded for n — 4 x lO -3 ^ 1 in the lower 
plate. 



5.4 Theory of asymptotic behavior for large g. 

We present a theory of the precipitation branch in the 
limit j > 1 to determine x p and n p . Assuming the 
branch (5.12) at the starting point, we confirm its ex- 
istence self-consistently. 

We first neglect the term —ksTgnp in Eq.(5.10) from 
gvonp <^ 1 and the term f(<j) a ) in Eq.(5.11) from <fi a = 1. 
In fact gvonp <C 1 in Fig. 9. We then obtain 

h = f{4>0) = -\f(4>p) + k B Tn a ]/(l - fa). (5.13) 

This determines the solute density n a in phase a as a 
function of <j)p in the form, 

n a S G(<i>p)/k B T, (5.14) 

where G(<j)) is a function of <f> defined as 

G(<f) = -f(<p)-(l-(f>)f(<t>) 

= -(fc B T/« o )[ln0 + X (l- ( / ) ) 2 ]. (5.15) 

From dG/d(f> = -(1 - 4>)f"{<j>) < and G(l) = 0, we 
have G(4>) > outside the coexistence curve, ensuring 
n a > in Eq.(5.14). 

In Eq.(5.10), we next use Eq.(3.3) for f'(<fi a ) to obtain 

Vo 1 k B T[- ln(l - <f> a ) - x - gv Q n a ] = f'(<f>p), (5.16) 

where the logarithmic term (oc ln(l — cf>)) balances with 
the solvation term (oc gn a ). Use of Eq.(5.15) gives 

l-4> a ^A p exp[-gv G(^)/k B T], (5.17) 

where the coefficient Ap is given by 

= exp[- x - vof'(<f>p)/k B T\, (5.18) 
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so Ap is of order unity. The factor exp[— gv G((j)p)/k B T] 
in Eq.(5.17) is very small for g >• 1, leading to 4> a — 1. 

Furthermore, from Eqs.(5.6) and (5.7), the volume 
fraction 7 Q of phase a is approximated as 



The above relation is rewritten as 



(5.19) 



G^) = k B Tn/(~f a 



-gA<j> 



k B Tn(l - 4>p) 



<f>-<f>p + (l- (j)p) exp[-5f(l 



(5.20) 



From the first to second line, we have used Eq.(5.7) and 
replaced A</> by 1 — (f>p. This equation determines </>p 
and 7 Q = (</> — (j>p)/(l — <j)p)- We recognize that G(<pp) 
increases up to Tne 9A ^ = Tne 9 ^ 1- ^ as 7 a — > or as 
(j)R 4>. In this limit it follows the marginal relation, 



G{#) = k B Tne 9{1 -^ {j a -> 0). 



(5.21) 



If n is fixed, this relation holds at \ — Xp so that 

Xp^I-ln^-wone^-^J/a-^) 2 , (5.22) 

where we use the second line of Eq.(5.15). Here n appears 
in the combination ne 9 ( 1-< ^ (3> n). On the other hand, if 
X is fixed, Eq.(5.21) holds at n = n p . Thus the minimum 
solute density n p is estimated as 



i p = e-^ x -®G{4>)/k B T, 



(5.23) 



velocity field v(r,t) being determined by the Stokes ap- 
proximation in Eq.(3.37). The cations and anions obey 



drij 
~dt 



V • (mv) 



D 



(5.24) 



where i = 1,2. The chemical potentials fii are defined 
in Eq.(3.11) and the ion diffusion constants are com- 
monly given by D. The space mesh size is a — v\ 
with d — 2. We measure time in units of to = a 2 Vo/L 
and set D = a 2 /to and 770 = 0.15/csT , to/woj where Lq is 
the kinetic coefficient for the composition and rjo is the 
shear viscosity. 

In Fig.ll, we show the time evolution of the droplet 
volume fraction, which is the fraction of the region <f> > 
0.6. In the early stage, the droplet number decreases 
in time with the evaporation and condensation mecha- 
nism. In the late stage, it changes very slowly tending 
to a constant. In our simulation without random noise, 
the droplets do not undergo Brownian motion and the 
droplet collision is suppressed. We also performed a sim- 
ulation for a neutral solute with g — 10.5 (not shown 
here), which exhibits almost the same phase separation 
behavior as in Fig.ll. 
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which is much decreased by the small factor e^ 9 ^ 1 ^^. oj 
In Fig. 10, the curves of x P and n p nearly coincide with _g 
the asymptotic formulas (5.22) and (5.23) in the range > 
4> < 0.35 for x p and in the range x < 2 for n p . They 
exhibit a minimum at (j> ~ e~ 9 /v ng for Xp and at ~ 
<7 _1 for n p . For larger <f> > 0.35, Xp nearly coincide with 
the coexistence curve, indicating disappearance of the 
precipitation branch. Notice that n p decreases to zero 
as 4> approaches the coexistence composition CX = 0.249 
at x = 2.2 (top) and 0.204 at x = 2 -3 (bottom), where 
phase separation occurs without solute. 

5.5 Simulation of spinodal decomposition for 
hydrophilic ions. In our theory^ we investigated 
solute-induced nucleation starting with homogeneous 
metastable states outside the spinodal curve (dotted 
line) in the left panels of Fig. 9. Here, we show two- 
dimensional numerical results of spinodal decomposition 
for hydrophilic ions with g\ = 12, 52 = 9, e 2 /eok B T = 
3a, and £\ = £0. At t = 0, we started with point (A) 
inside the spinodal curve in the left lower panel in Fig. 9 
using the common values of the static parameters given 
by x = 1-92, 4> = 0.58, and n = (m + n 2 ) = 4 x 10~ 3 /«o- 
In addition, we set C — 2k B Ta 2 /v . On a 256 x 256 
lattice under the periodic boundary condition, we in- 
tegrated Eq.(3.36) for the composition <j>(r,i) with the 




0.15 



0.05 




t = 24002St 
1 



t = 108026tt, ( - 180026t 



10 4 */£o 



0.5 



FIG. 11: Time evolution of droplet volume fraction (that of 
the region cj> > 0.6) and composition snapshots of precipitated 
droplets at three times. They are induced by hydrophilic ions 
with g\ = 12 and gi = 9. The initial state was at point (A) 



in Fig.9 with x = 1-92, cj> = 0.58, and n = 4 x 10" d 
droplet volume fraction nearly tends to a constant. 



The 



VI. THEORY OF POLYELECTROLYTES 

6.1 Weakly ionized polyelectrolytes. In this sec- 
tion, we consider weakly charged polymers in a theta 
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or poor, one-component water-like solvent in the semidi- 
lute case <p > A -1 / 2 . Following the literature of poly- 
mer physics^, we use and N to represent the poly- 
mer volume fraction and the polymerization index. Here 
charged particles interact differently between uncharged 
monomers and solvent molecules. The selective solvation 
should become more complicated for mixture solvents, as 
discuused in Sec.l. 

To ensure flexibility of the chains, we assume that the 
fraction of charged monomers on the chains, denoted by 
/i on , is small or / ion <C 1. From the scaling theory^, the 
polymers consist of blobs with monomer number g b = 
4>~ 2 with length ^ = ag\J 2 — ac/)^ 1 . The electrostatic 
energy within a blob is estimated as 

e b = k B T(f ion g) 2 £ B /t; b = k B Tf 2 J B /^ 3 a. (6.1) 

where l B is the Bjerrum length. The blobs are not much 
deformed under the weak charge condition ej, < k B T, 
which is rewritten as 



•2) 



6.2 Ginzburg- Landau theory. The number of the 
ionizable monomers (with charge — e) on a chain is writ- 
ten as vmN with vm < 1- Then the degree of ionization 
(or dissociation) is C = /km/^M and the number density 
of the ionized monomers is 



n p = v Vion^ = v o 1 vm(,4>, 
The charge density is expressed as 

p = e ^ Z i n i ~ en p- 

i=c, 1,2 



(6.3) 
(6.4) 



Here i = c represents the counterfoils, i = 1 the added 
cations, and i = 2 the added anions. The required rela- 
tion /; on <C 1 becomes vmC <C 1 (which is satisfied for 
any ( if v M < 1). 

We set up the free energy F accounting for the molec- 
ular interactions and the ionization equilibrium^!. Then 
F assumes the standard form (3.1), wh ere t he coefficient 
of the gradient free energy is written a d 5 * 36 ! 



C(</>) = k B T/l2a(j>{l 



■5) 



1/3 

in terms of the molecular length a = v and (j). The 
/tot consists of four parts as 

/tot = f{<t>) + k B T ruMniXf) - 1 + g i( f>\ 

i=c,l,2 

+k B T(A a + g p 4>)n p + f dis . (6.6) 
The first term / is of the Flory-Huggins forrrpMI, 
v f 



k B T N 



In0+(l-0)ln(l-0)+X<^(l-0). (6-7) 



The coupling terms (oc gi,g p ) arise from the molecular 
interactions among the charged particles (the ions and 



the charged monomers) and the uncharged particles (the 
solvent particles and the uncharged monomers), while 
k B T A is the dissociation free energy in the dilute limit 
of polymers (<f> — > 0). The last term in / tot arises from 
the dissociation entropy on chainJ*nHHl 



_vo_ 
k B T 



/dis — Vm<, 



ClnC+(l-C)Ml-C) 



6.3 Dissociation equilibrium. If F is minimized 
with respect to (, it follows the equation of ionization 
equilibrium or the mass action law, 



C C/(1 ^ C) = K (</>), 



(6.9) 



where n c is the counter ion density and K(<f>) is the dis- 
sociation constant of the form, 



K{4>) = v Q 1 exp[-A - (g p + g c )(f>}. 



(6.10) 



We may interpret k B T[A + (g p + g c )(f>] as the 
composition-dependent dissociation free energy. With in- 
creasing the polymer volume fraction (j), the dissociation 
decreases for positive g p + g c and increases for negative 
9 P + 9c- If g P + ffc ^ 1, K{4>) much decreases even for a 
small increase of 4>. Here K((f>) has the meaning of the 
crossover counterion density since £ is expressed as 



C=l/[l+n c /A'(0)] 



(6.11) 



which decreases appreciably for n c > K(cj>). 

In particular, if there is no charge density and no salt 
(rtp = n c and ni = ria = 0), n c satisfies the quadratic 
equation n c (n c + K) = v^ 1 i>m4>K, which is solved to give 



C = v a n c /v M (t) = 2/(v / Q(</>) 
Here it is convenient to introduce 



1 + 1 



3.12) 



(6.13) 



We find C < 1 and n c = {vm&K/vq) 1 / 2 for Q > 1, while 
£ — > 1 for Q <C 1. The relation (6.12) holds approxi- 
mately for small charge densities without salt. 

6.4 Structure factor. As in Sec. 3, it is straightfor- 
ward to calculate the structure factor S(q) for the fluc- 
tuations of (f> on the basis of /tot in Eq.(6.6). As a func- 
tion of the wave number q, it takes the same functional 
form as in Eq.(3.14), while the coefficients in the poly- 
electrolyte case are much more complicated than those in 
the electrolyte case. That is, the shift — (gi + g-2) 2 no/2 in 
Eq.(3.14) is replaced by its counterpart Ar dependent on 
ni and n p (for which see our paper^l). In the following 
expressions (Eqs.(6.14)-(6.17)), m, n p , and (j) represent 
the average quantities. The Debye wave number of poly- 
electrolytes is given bj®l 



k 2 = Att£ b 



(l-()n p + J2 Z ? n i 



(6.14) 
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which contains the contribution from the (monovalent) 
ionized monomers (oc n p ). The asymmetry parameter 7 P 
in Eq.(3.14) is of the form, 

7p = (47r£ B k B T/C)^ 2 A/ K 2 , (6.15) 
where C is given by Eq.(6.5) and 

A = n p /cj>- (1 - C)g p n p + '^2,Z i g i n i . (6.16) 

i 

In this definition, 7 P can be negative depending on the 
terms in A. Mesophase formation can appear for |7 p | > 1 
with increasing \- 

The parameter 7 P is determined by the ratios among 
the charge densities and is nonvanishing even in the dilute 
limit of the charge densities. In particular, if n c = n p and 
n i = n 2 in the monovalent case, 7 p is simplified as 



7 P 



cp- 1 -{l-Qg p + g c + { gi -g 2 )R 
(An£ B C/k B T) 1 / 2 (2-C + 2R) ' 



(6.17) 



where the counterions and the added cations arc differ- 
ent. The R = ni/n p is the ratio between the salt den- 
sity and the that of ionized monomers and Eq.(6.5) gives 
(47r£ B C/k B Ty/ 2 = [*e B /3a<f>(l - dp)] 1 ' 2 . 

Some consequences follow from Eq.(6.17). (i) With 
enriching a salt we eventually have R \gi\; then, the 
above formula tends to Eq.(3.15), which is applicable for 
neutral polymer solutions (and low-molecular-weight bi- 
nary mixtures for N = 1) with salt, (ii) Without the 
solvation or for gi = 0, the above S (q) tend s to the previ- 
ous expressions for polyelectrolytes^ 1 931941 , where 7 P de- 
creases with increasing R. In accord with this, Braun 
et alPSl observed a mesophase at low salt contents and 
macrophase separation at high salt contents, (iii) In our 
theory, neutral polymers in a polar solvent can exhibit a 
mesophase for large \g± — g 2 \- 

Hakim et a/! 96 * 97 ! found a broad peak at an intermedi- 
ate wave number in the scattering amplitude in (neutral) 
polyethylene-oxide (PEO) in methanol and in acetoni- 
trile by adding a small amount of salt KI. They ascribed 
the origin of the peak to binding of K + to PEO chains. 
Here more experiments are informative. An experiment 
by Sadakane et alW^ suggests that use of an antagonistic 
salt would yield mesophases more easily. 

6.5 Interface profiles without salt. We suppose 
coexistence of two salt-free phases (m = n 2 = 0), sepa- 
rated by a planar interface. Even without salt, the inter- 
face profiles are extremely varied, sensitively depending 
on the molecular interaction parameters, Ao, g p , and g c . 
If a salt is added, they furthermore depend on g x , g 2 , and 
the salt amount. The quantities with the subscript a (/3) 
denote the bulk values in the polymer-rich (solvent-rich) 
phase attained as z — > — oo (as z — ¥ oo). The ratio of the 
bulk counterion densities is given by 



exp 



g c A(/) 



eA$ 



k R T 



(6.18) 



The Galvani potential difference A$ = $ Q — is ex- 
pressed in terms of Q(<j>) in Eq.(6.13) as 



eA$ 



g p A4> + In 



y/QWfi 



1 - 1 



+ 1-1 



(6.19) 



If Q(<t> a ) > 1 and Q{ct>f)) > 1 (or Cq < 1 and C/3 < 1), 
we obtain eA$/k B T = (g p - g c )Acf)/2 + )xi{(j>p/(j> a ). 
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FIG. 12: Interface profiles in the salt-free case for (a) Ao = 5, 
g p = 1, and g c = 4 (top) and for (b) Ao = 8, g p = 2, and 
g c — —6 (bottom). Polymer volume fraction (p(z), normalized 
potential 4>(z)/$m, and degree of ionization £(z) (left), and 
normalized charge densities von c (z) and von p (z) (right). The 
other parameters are common as x = 1> N — 20, vm = 0.1, 
£i = — 0.9eo, and Ib = 8a/n. Here <&(z) is measured from 
its minimum, and $m(= 2.67fcsT/e in (a) and 6.86fc_BT/e in 
(b)) is the difference of its maximum and minimum. [From: 
A. Onuki and R. Okamoto, J. Phys. Chem. B, 113, 3988 
(2009).] 



We give numerical results of one-dimensional profiles in 
equilibrium, where we set x = 1, iV = 20, E\ — — 0.9eo> 
and l B = e 2 /e T = 8a/7r. The dielectric constant of 
the solvent is 10 times larger than that of the polymer. 
The space will be measured in units of the molecular size 

1/3 

a = v . In Fig. 12, we show salt-free interface profiles 
for (a) A = 5, g p = 1, and g c — 4 and (b) A = 8, 
g p = 2, and g c = — 6. In the a and /? regions, the degree 
of ionization £ is 0.071 and 0.65 in (a) and is 0.24 and 
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FIG. 13: Periodic profiles in a salt-free mesophase with vu = 
0.5, A = 5, g p = 1, and g c = 4. Top: (j>(z). $(z)/$ M 
with $m = 0.86T/e, and (,(z). Bottom: von c (z) and von p (z). 
[From: A. Onuki and R. Okamoto, J. Phys. Chem. B, 113, 
3988 (2009).] 



0.51 in (b), respectively. The normalized potential drop 
e($ a - § p )/k B T is -2.67 in (a) and 0.099 in (b). In- 
terestingly, in (b), Q(z) exhibits a deep minimum at the 
interface position. We can see appearance of the charge 
density n c — n p around the interface, resulting in an elec- 
tric double layer. The counterion density is shifted to 
the (3 region in (a) because of positive g c and to the a 
region in (b) because of negative g c . The parameter 7 P 
in Eq.(6.16) is 0.75 in (a) and 0.20 in (b) in the a region, 
ensuring the stability of the a region. 

The surface tension a is again expressed as in Eq.(3.26) 
with the negative electrostatic contribution. It is calcu- 
lated as a = 0M75k B T/a 2 in (a) and as 0.0556fc B T/a 2 
in (b), while we obtain a = 0.050fcsT/a 2 without ions at 
the same X = 1. In (a) a is largely decreased because the 
electrostatic term a c in Eq.(3.27) is increased due to the 
formation of a large electric double layer. In (b), on the 
contrary, it is increased by 10% due to depletion of the 
charged particles from the interface^. 

We mention calculations of the interface profiles in 
weakly charged polyelect rolyte s in a poor solvent using 
self-consistent field theorjJ2§l22l. In these papers, however, 
the solvation interaction was neglected. 

6.6 Periodic states without salt. With varying 
the temperature (or the average composition (4>), the 
amount of salt, there can emerge a number of mesophases 



sensitively depending on the various molecular parame- 
ters (<?j, Aq, and I'm)- In Fig. 13, we show an example 
of a one-dimensional periodic state without salt. Here 
vm is set equal to 0.5 and the charge densities are much 
increased. In this case, the degree of segregation and the 
charge heterogeneities are much milder than in the cases 
in Fig. 12. 



VII. SUMMARY AND REMARKS 

In this review, we have tried to demonstrate the crucial 
role of the selective solvation of a solute in phase tran- 
sitions of various soft materials. We have used coarse- 
grained approaches to investigate mesoscopic solvation 
effects. Selective solvation should be relevant in under- 
standing a wide range of mysterious phenomena in wa- 
ter. Particularly remarkable in polar binary mixtures are 
mesophase formation induced by an antagonistic salt and 
precipitation induced by a one-sided solute (a salt com- 
posed of hydrophilic cations and anions and a neutral hy- 
drophobic solute). Regarding the first problem, our the- 
ory is still insufficient and cannot well explain the compli- 
cated phase behavior disclosed by the experimentiP^itll] 
To treat the second problem, we have started with the 
free energy density /tot in Eq.(5.1), which looks rather 
obvious but yields highly nontrivial results for large g. 
Systematic experiments are now possible. In particu- 
lar, this precipitation takes place on colloid surfaces as a 
prewetting phase transition near the precipitation curve 
X = X P as in Fig.lcP. 

Though still preliminary, we have also treated an ionic 
surfactant system, where added in water-oil are cationic 
surfactant, anionic counterions, and ions from a salt. In 
this case, we have introduced the amphiphilic interaction 
as well as the solvation interaction to study the inter face 
adsorption. For ionic surfactants, the Gibbs formulaP^ 
for the surface tension is insufficient, because it neglects 
the electrostatic interaction. 

In polyelectrolytes, the charge distributions are ex- 
tremely complex around interfaces and in mesophases, 
sensitively depending on the molecular interaction and 
the dissociation process. Our continuum theory takes 
into account these effects in the simplest manner, though 
our results are still fragmentary. Salt effects in polyelec- 
trolytes should also be further studied, on which some 
discussions can be found in our previous papep^Sl In the 
future, we should examine phase separation processes in 
polyelectrolytes, where the composition, the ion densi- 
ties, and the degree of ionization are highly inhomoge- 
neous. In experiments, large scale heteroge neities have 
been observed to be pinned in space and tim d 100 * 1 ^, giv- 
ing rise to enhanced scattering at small wave numbers. 

As discussed in Sec.l, there can be phase separation 
induced by selective hydrogen bonding. In particular, 
the effect of moisture uptake is dramatic in PS-PVMEp^, 
where scattering experiments controlling the water con- 
tent are desirable. To investigate such polymer blends 
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theoretically, we may use the form in Eq.(5.1) with n be- 
ing the water density and /(</>) being the Flory-Huggins 
free energy for polymer blends^. Similar problems should 
also be encountered in block polymer systems containing 
ions or water. It is also known that blends of block poly- 
mer and homopolymer exhibit complicate d ph ase behav- 
ior for different interaction parameters Xij^- 

We mention two interesting effects not discussed in this 
review. First, there can be an intriguing interplay be- 
tween the solvation and the hydrogen bonding in phase 
separation. For example, in some aqueous mixtures, even 
if they are miscible at all T at atmosphere pressure with- 
out salt, addition of a small amount of a hydrophilic salt 
gives rise to reentrant phase separation behavior^EH 
On the other hand, Sadakane et al. observed a shrinkage 
of a closed-loop coexistence curve by adding an atago- 
nistic salt or an ionic surfactant!^. Second, molecular 
polarization of polar molecules or ions can give rise to 
a surface potential difference on the molecular scal e at 
an interface. See such an example for wat er-hexan d 1 3 1 
This effect is particularly noteworthy for hydronium ions 
in acid solution^. 

We will report on the wetting transition on charged 
walls, rods, and colloids and the solvation-induced col- 
loid interaction. These effects are much influenced by 
the ion-induced precipitation discussed in Sec. 5. In these 
problems, first-order prewetting transitions occur from 
weak-to strong ionization and adsorption, as discussed in 
our paper on charged rodiP^. We will also report that a 
small amount of a hydrophobic solute can produce small 
bubbles in water even outside the coexistence curve, on 
which there have been a large number of experiments. 
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Appendix A: Statistical theory of selective 
solvation at small water composition 

We present a simple statistical theory of binding of 
polar molecules to hydrophilic ions due to the ion-dipole 
interaction in a water-oil mixture^ when the water vol- 
ume fraction <fi is small. We assume no macroscopic inho- 
mogeneity and do not treat the large-scale electrostatic 
interaction. Similar arguments we re giv en for hydrogen 
boning between water and polyme r 54 * 55 ! . 

Our system has a volume V and contains N w water 
molecules. Using the water molecular volume vq, we have 



4> = N w v /V = N w /N , 



(Al) 



where No = V/vq. We fix A w or <j> in the following. The 
total ion numbers are denoted by Nu — Vrii , where i = 1 
for the cations and i = 2 for the anions with rii being the 



average densities. The ionic volumes are assumed to be 
small and their volume fractions are neglected. Then the 
oil volume fraction is given by 1 — <j>. Each solvation shell 
consists of v water molecules with v — 1, • • • , Zi, where 
Zi is the maximum water number in a shell. 

Let the number of the J/-clusters composed of v water 
molecules around an ion be 7i„-/V w . The total number of 
the solvated ions is then "fiNu with 



7i 



,liu < 1, 



(A2) 



where 1 < v < Zi. The number of the bound water 
molecules in the ^-clusters is vjwNn. The fraction of 
the unbound water molecules 4>f satisfies 



(A3) 



We construct the free energy of the total system F tot 
for each given set of ji U . In terms of the oil density n ji = 
Vq 1 (1 ~ <j)), the unbound water density n w f = Vq 1 ^, the 
unbound ion densities na = rii(l — ji), and the cluster 
densities m v = n^w, we obtain 



F, 



Vk B T 



n oi i[ln(n ol iAo 11 ) - 1] + n w{ [ln(n wf A 3 



"if 



ln(n if Af 



\ " 



//, +> n iu [ln(n iv \f u ) - w° v ] 

(A4) 



where A n, A w , Ai, and \i V are the thermal de Broglie 
wavelengths, ksTw^ are the "bare" binding free ener- 
gies, and x is the interaction parameter between the un- 
bound water and the oil. We asssume short-range inter- 
actions among the clusters and the oil characterized by 
the interaction parameters x%v to obtain the last term. 
At small <p, the interactions among the clusters and the 
unbound water are neglected. That is, we neglect the 
contributions of order (f> 2 . We then calculate the solva- 
tion contribution F so \ = F tot — Fq, where Fq is the free 
energy without binding (7^ = 0). Some calculations give 



F so \ 
N k B T 



= &(ln&-l)+5^TH(l-7i)ln(l-7i) 



+ y"l ? yfo/( lll 7»/ - Wiv) ~ </>(ln<£- 1), (A5) 

where ksTw^ are the "renormalized" binding free ener- 
gies written as 



w° v +31n(AiA^/Ai„) +vx~ X%v 



(A6) 



The fractions 7^ are determined by minimization of _F so i 
with respect to 7^ under Eqs.(A2) and (A3) as 



7j = i-i/(i+5>r^ 



(A7) 
(A8) 
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Substitution of Eqs.(A7) and (A8) into Eq.(A5) yields 
F so \ 



k B TN 



= <t>\n^ + <f> - fa + ^«oniln(l-7i). (A9) 



First, we assume the dilute limit of ions N] <C N w , 
where we have cj> — <C <p- In the right hand side 
of Eq.(A9), the sum of the first three terms becomes 
— (4> — 4>f) 2 /2<fi and is negligible. We write F so i as the 
sum V^2 i riinl ol ((j>) and use Eq.(A8) to obtain the solva- 
tion chemical potentials of the form, 

/4iW = -k B Tln(l + J2re w ^). (A10) 



Let the maximum of k B Twi^/v (per molecule for various 
v) be €hi for each i (see Eq.(2.3)). Then we obtain the 
expression (2.2) for the crossover volume fraction <^* ol . 
For (j> > 4>l ol , 7j approaches unity. 

Second, we consider the dilute limit of water, <fi <C 
4>l ol (i = 1,2), where the cluster fractions i/< are small 
and the dimers with v = 1 are dominant as indicated in 
the experiment!^. Neglecting the contributions from the 
clusters with v > 2, we obtain 



(All) 



(A12) 



Ti^&e^ 1 , f = 0/(1 + 5'), 

where S is the parameter defined as 

S = v a (me W11 +n 2 e W21 ). 

The solvation free energy behaves as 

F sol /N o k B T = -0\n(l + S). (A13) 

For S < 1, we find MLlW - -k B T<j)e Wl '. However, 
if S > 1, the solvation chemical potential are not well 
defined. 

Appendix B: Ions at liquid-liquid interface 

In electrochemistry, attention has been paid to the ion 
distribution and the el ectric potential difference across 
a liquid-liquid interfac e! 75 * 76 !. (In the vicinity of an air- 
water int erface , virtually no ions are present in the bulk 
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air regio n * I) Let us suppose two species of ions (i = 
1,2) with charges Z\e and Z 2 e {Z\ > 0, Z 2 < 0). At 
low ion densities, the total ion chemical potentials /ij in 
a mixture solvent are expressed as 



where is the thermal de Broglie length (but is an ir- 
relevant constant in the isothermal condition) and $ is 
the local electric potential. This quantity is a constant in 
equilibrium. For neutral hydrophobic particles the elec- 
trostatic term is nonexistent, so we have Eq.(2.5). 

We consider a liquid-liquid interface between a polar 
(water-rich) phase a and a less polar (oil-rich) phase ft 
with bulk compositions <fi a and 4>p with <p a > cbp. The 
bulk ion densities far from the interface are written as 
n ia in phase a and n^p in phase ft. From the charge 
neutrality condition in the bulk regions, we require 



Zinia + Z 2 n 2a = 0, ZiTiip + Z 2 n 2 fj = 0. 



(B2) 



The potential <£> tends to constants $ a and $^ in the 
bulk two phases, yielding a Galvani potential difference, 
A$ = $ a — Here $ approaches its limits on the 
scale of the Debye screening lengths, k^ 1 and k^ 1 , away 
from the interface, so we assume that the system extends 
longer than k" 1 in phase a and k^ 1 in phase ft. Here we 
neglect molecular polarization of solvent molecules and 
surfactant molecules at an interface (see comments in the 
summary section). 

The solvation chemical potentials Msoi(^) a -^ so take dif- 
ferent values in the two phases due to their composition 
dependence. So we define the differences A/x^ as in 
Eq.(2.4). The continuity of /Xj across the interface gives 

k B T\n(n la /n lP ) + Z,eA$ - A^ = 0, (B3) 

where i — 1,2. From Eqs.(B2) and (B3), the Galvani 
potential difference A$ is expressed a a 23 * 76 ! 

A$ = [A^-A^]/ e (Z 1 + |Z 2 |). (B4) 

Similar potential differences also appear at liquid-solid 
interfaces (electrodes)^. The ion densities in the bulk 
two phases (in the dilute limit) are simply related by 



nip _ "2/3 
nia n 2a 



exp 



\Z 2 \)k B T 



(B5) 



IH = fc B Tln(n I A- ) + Z t e9 + /4i(<W, 



(Bl) 



However, if three ion species are present, the ion parti- 
tioning between two phases is much more complicated^. 
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